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Abstract 

Existence of amplitude independent frequencies of oscillation is an unusual 
property for a nonlinear oscillator. We find that a class of N coupled nonlin- 
ear Lienard type oscillators exhibit this interesting property. We show that 
a specific subset can be explicitly solved from which we demonstrate the ex- 
istence of periodic and quasiperiodic solutions. Another set of iV-coupled 
nonlinear oscillators, possessing the amplitude independent nature of fre- 
quencies, is almost integrable in the sense that the system can be reduced to 
a single nonautonomous first order scalar differential equation which can be 
easily integrated numerically . 

Keywords: Nonlinear oscillators, coupled ordinary differential equations, 
complete integrability, isochronous systems 

1. Introduction 

Let us consider the following nonlinear oscillator described by the so called 
modified Emden equation with linear forcing term [lj], 

x + 3xx + x 3 + u 2 x = 0. (1) 

Here u is a parameter. Equation (1) can be considered as the cubic anhar- 
monic oscillator with additional position dependent damping type nonlinear 
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force 3xx. This type of equation has been well studied in the literature. For 
example, Eq. (1) with u = arises in a wide range of physical problems: 
it occurs in the study of equilibrium configurations of a spherical gas cloud 
acting under the mutual attraction of its molecules and subject to the laws of 
thermodynamics 0, l^] and in the modelling of the fusion of pellets 0] ■ It also 
governs spherically symmetric expansion or collapse of a relativistically grav- 
itating mass [5] . This equation can also be thought of as a one-dimensional 
analog of the boson gauge-theory equations 0, 0]. 

Equation flTJ has been shown to posses an unusual property which is not 
a general characteristic of a nonlinear equation: The frequency of oscillation 
of the oscillator is independent of the amplitude similar to that of a linear 
harmonic oscillator An oscillator which possesses this property is also 
known as an isochronous oscillator [8J. For a detailed study about isochronous 
orbits and isochronous oscillators one may refer to Refs. 0, [9| . Equation ([1]) 
admits the following nonsingular, periodic solution: 



A sin (ut + S) 
[l-£cos(ut+6)) 



x(t) = - "T. _,Z ^ . 0<A<u. (2) 



Here A and 5 are arbitrary constants, expressible in terms of the two integrals 
of motion or integration constants obtained by solving (JTJ (for details see ref. 

Note that the angular frequency of oscillation u continues to be the same 
as that of the linear oscillation. From this solution it is obvious that for < 
A < u, equation ([T]) exhibits the property of amplitude independence of the 
frequency of oscillation. One can starightforwardly write down the solution 
of the initial value problem from the general solution (j2J). For example, for 
the initial condition x(0) = B = , A 2 , x(0) = 0, from ([2]) we have the 

solution as 



Bu sin 



ut + cos 1 ^ 



B 



*(*) = , — r — 7 - ^ - (3) 

V B 2 + u 2 — B cos ut + cos^ 1 ( - 



Note that B is the amplitude of oscillation. Figure [1] shows the periodic 
oscillations admitted by Eq. (CQ) for three different sets of initial conditions 
x(0) and x(0) with u = 1.0 in terms of three different colours. We note here 
that the frequency of the oscillations is independent of the initial conditions 
as in the case of the linear harmonic oscillator. 



2 




10 20 30 40 -0.8 -0.4 0.4 0.8 

t X 

Figure 1: (Color online) (a) Time series plot of Eq. ([I]) exhibiting periodic oscillation 
for three different initial conditions (three different colours) and to — 1.0 (b) Phase space 
portrait of Eq. (fT]) 



One can trace the origin of this property of equation (TT]) to the fact that 
it can be transformed to the linear harmonic oscillator equation, 

U + uj 2 U = 0, (4) 

through a nonlocal transformation, 

U = x{t)eti x(t ' )dt '. (5) 

The solution fl2]) can be obtained (see below, equation (fT2]) ) from the solution 
of fll]), U = Asm(ut + 5), where A and 5 are arbitrary constants and the fre- 
quency, co, is independent of the amplitude. Such a linearization property is 
one of the fascinating features associated with a class of nonlinear equations 
exhibiting large number of symmetries and extensive search for such lineariz- 
ing transformations is being made in the recent literature [H, [lo|, 11, 12]. In 



fact, there exists a class of nonlinear oscillators which are connected to the 
linear oscillator equation (j3]) through the following nonlocal transformation 



12 



U = x(t)eX f W»'*, (6) 
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where f(x(t)) is an arbitrary function of x(t). Now substituting (jH]) into (jl]) 
we get a nonlinear ordinary differential equation (ODE) of the form 

x + (2f + xf')x + (f 2 + u 2 )x = 0, (7) 

where prime denotes differentiation with respect to x. Equation (JJJ) is a 
special case of the well known Lienard equation (LE) [l3[ 

x + u(x)x + v(x) = 0. (8) 

One can also consider a more general nonlocal transformation of the form 

U = g(x(t))ef , (9) 

and substituting this in (Jl]) we get 

i + il^ + il i+2fi + ^i + ti = , ( gl = f). (10) 

9 9 9 9 V dx J 

We find the above equation reduces to a Lienard type equation only for the 
choice g(x) = x. Interestingly for / = x~ 2 , equation (J7|) becomes the well 
known isotonic oscillator 0] equation, 

X + U 2 X+ 4r = 0. (11) 
X 6 

The solution of the nonlinear equation ([7]) is obtained by using the identity 

^ = £ + /(*(*)). (12) 
U x 

Since U = A sin (ut + 5) , where A and 5 are integration constants, is the 
solution of the linear harmonic oscillator (jl]), equation (fl2]l can be rewritten 
as the first order nonlinear differential equation of form 

OJX 

x + f(x(t))x 7 ft = 0. (13) 

J K y " tan(wt + 5) v ; 

Now one can get the solution of ((Tj) by solving (jTBl) . In particular, for the 
specific case f = x g equation (fl3"|) becomes a Bernoulli equation of the form 




Figure 2: (a) Time series plot of Eq. (IT51) for q — 3 exhibiting periodic oscillations with 
the initial condition x(0) = 3 and i(0) = for u = 1 (b) Phase space portrait of Eq. (jl~5l) 



The corresponding ODE (J7J) becomes 

x + (q + 2)x q x + x 2q+l + u 2 x = 0, (15) 

and equation ([I]) is the special case corresponding to q — 1. 
Upon integrating (Tl4|) we get the periodic solution of (|T5|) as 

*W = ^ .-, (16) 



[I - cos(t)(sin 2m (t) + J2T=o C k sin 2m - 2fe - 2 (F)] 



(2ni+l) 



771- 



wheret = wt + 5, |^ =0 C fe |</-l, g = 2m + l, C fe = (2 fc+1 m(m-l) 
k))/((2m — l)(2m — 3) . . . (2m — 2k — 1)), J and 5 are arbitrary constants. 
Here m is a non-negative integer and u is the angular frequency. One can 
note that solution (Tl6|) is also isochronous. This has indeed been reported 
recently by Iacono and Russo 1(| using a different procedure. In figure [2] we 



show the periodicity of the solution for the case q = 3 and with the initial 
conditions x(0) = 3 and x(0) = 0. We additionally remark here that the 
case q = 2m, m = 1, 2, . . . , of equation ffl5|) is also exactly solvable but the 
solutions are of damped oscillatory type as will be proved later in this article 
(Sec. [2]) even for coupled systems. 

In this paper we will show that this unusual (amplitude independent 
frequency) property possessed by the class of equations (J7|) is not limited 
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to scalar LEs alone but is also possessed by a class of N coupled nonlinear 
oscillator equations. In order to demonstrate the existence of periodic and 
quasiperiodic solutions for such systems, we extend the above mentioned 
procedure to N coupled LEs. We also derive the integrals of motion and the 
general solution to the system of N coupled LEs. We also point out that 
there exists another class of N coupled nonlinear oscillators which is almost 
integrable in the sense that it admits (2N — 1) independent integrals and 
reduces to a single first order nonautonomous nonlinear differential equation 
that can be solved numerically. The frequency of oscillation of its periodic 
solutions is again independent of amplitude. 

In the following section we show that a system of N coupled LEs can 
be explicitly integrated to find periodic and quasiperiodic solutions whose 
frequencies are identical to a system of uncoupled linear oscillators. In Sec. 
|3]we deduce the underlying 2N integrals of motion for the system and analyze 
their structure. In Sec. H]we numerically solve the system of iV coupled LEs 
for which the explicit general solution is not known, but (2N — 1) integrals 
are known. We finally summarize our results in Sec. 

2. Coupled nonlinear oscillators of Lienard type and integrability 

Generalizing the above results of the scalar system (j7|) to a system of 
coupled nonlinear oscillators of Lienard type, we relate them to a system of 
uncoupled iV-dimensional linear harmonic oscillators. For this purpose, we 
consider the N- dimensional anisotropic harmonic oscillator equation of the 
form 

Ui + ufUi = 0, z= 1,2,..., TV, (17) 

where the freqencies u i: i = 1,2, . . . , N are in general different. Let us now 
introduce the nonlocal transformation 

Ui = XlC i ^^-x N )dt = x . e Ifi(x)d tj i = 1,2, . . . ,N, (18) 

where fi(x) = fi{x\, X2, ■ ■ ■ xn) are arbitrary functions of the variables. Sub- 
sisting ( JIB]) into (JIT]) we get a system of N coupled Lienard type second 
order nonlinear oscillator equations of the form 

N 

Xi ~\~ ^ fi {x^XiXj ~\- fi{x^Xi 

+f?(x) Xi + u? Xi = 0, % = 1, 2, . . . , N, (19) 
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where = dfi/dxj, i,j — 1, 2,...,N. In the following we shall demon- 
strate the existence of periodic and quasiperiodic solutions with amplitude 
independent frequencies, u iy % = 1,2, ... ,N (isochronous property) for the 
above mentioned nonlinear oscillator equation with specific forms of /Vs. 
To obtain the desired results, we make use of the following identities, 

^ = - + m, i = l,2,...,N, (20) 

derived from equation (|T8|) . Obviously the solution of the system of second 
order linear ODEs (jT7|) is 

Ui = At sm(uit + 5i), (21) 

where A^s and S^s, i = 1, 2, . . . , N, are 2N integration constants. Now 
substituting the solution ( l2Tj) in equation (|20|) we find the following system 
of coupled first order ODEs to represent (fT9l) . 

^ = , rTXT^ - ^ = 1, 2, . . . , iV. (22) 

It may be noted that N integration constants, 5iS, i = 1, 2, . . . , N, of the iV 
coupled second order ODEs ( fT9l) appear explicitly in ( !22l . In order to find 
the general solution of ( |T9|) we need A" more integration constants which are 
to be obtained by integrating ( 1221) . For general forms of fi(x) in ( 1T9]) this 
cannot be done. However for the special choice, 

fk = /jv(^) = /(^) = f(x 1 ,x 2 , ...x N ), k = 1, 2, . . . , N - 1, (23) 

equation (1221) can be integrated. In order to perform this integration we 
multiply each equation of the system ([22]) by xn for i = 1,2, . . . , N — 1 and 
the last equation by Xi and subtract them to get 



(■ 



x k x N - x N x k = [ - — i—^ — ^ - - — t-^—T\ ) x k%N, k = 1,2,..., N.(24) 
\ta.n{oj kt + Ok) ta,n{ujNt + On) 



(25) 



Dividing Eq. (J24[) throughout by XkXN, we get 

Xk Xn\ i ^k % 



Xk xn J \tan(cujtt + 5k) ta,n(u>Nt + 5 



N, 
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Integrating the above equation, we get 



logx k — logXjv = logsin(o; fc t + 5) — logsin(u;jv£ + 8) + log/ fc . (26) 
Rewriting the above equation, we get 
x k sin(u k t + 5 k ) 



xn sin(o;jvt + <5/v) 



h k (t), k = l,2,...,N-l (27) 



where Ik's, k — 1, 2, . . . , N — 1, are N — 1 integration constants. We are also 
left with a first order ODE of the form 

±N = I — HVt x \ ~ f(Ht),x N )x N , (28) 

where hit) = {hi(t), h,2(t), . . . , /tAr-i(t)} and hk(t)'s are given by equation 
( 1271) . Now the problem of solving the set of N coupled autonomous second or- 
der ODEs ( fl9l) is reduced to the problem of solving a single non-autonomous 
first order ODE (|28j) . Therefore one can get the general solution for equation 
(fl9|) for the case (|23|) by solving equation ( j28|) . 

Equation ([281 again j^annot be in general solved explicitly for arbitrary 



form of the function / |14j. Hence, in order to solve (1281) we assume that the 
function / has a symmetry f(oex) = a Q f(x), where ax = {axi, ax<i, ■ ■ ■ ax^}, 
a and q are arbitrary parameters. This implies that / is a homogeneous 
polynomial and we assume the following form of /, 

/ = I>(*K, (29) 
i=i 

where 7z(t)'s are arbitrary functions of t. However, in the present case we 
assume 7;'s to be constants only for simplicity. Even when they are functions 
of t the following integration procedure holds good. With the above choice 
of /, equation (Tl9l) reduces to the system of coupled nonlinear oscillator 
equations, 

N N 
Xi + 2 22 lltfXi + (1 X] ^iXitfxi 

1=1 1=1 

N 

+ EW) 2 ^ + ^ = 0, i = 1,2,..., N. (30) 
i=i 
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The solution of this system of coupled nonlinear oscillators can be obtained 
by solving the following first order nonlinear differential equation obtained 
by substituting (J2SD into (J2E5 along with ([27]): 



x N = u N cot{u N t + b N )x N - x N ) 7,— rT~TT ' ( 31 ) 

^ sm^w^t + 5jv) 

where 77 = Iiji, I = 1, 2, . . . , N — 1, and = 77V • The above equation is of 
the first order Bernoulli equation type 13j, ll4(, namely 

/77; 

— = P(u)v + Q(u)v n , (32) 
aw 

where -P(w) and Q(u) are arbitrary functions of the independent variable u. 
With the substitution y{t) = XAr(t) 9 in f l3T]) we get the following first order 
linear inhomogenous ODE 

N 



dy / e- \ - sm g (uit + 5;) 
— = -qu N cot(u N t + 5 N )y + q> 7« — "777 7 , x s 

(it -f— ' Sin (c^Tyt — )— TV y 



The general solution of 032]) is obviously 

y(f) = + 6n) , 04) 

where In is the integration constant. Rewriting the above in terms of x^, 
we get 

Mt) = Si \ iUNt + ^ , (35) 

where In is the 2N integration constant which we are looking for. 

The integral appearing in the denominator of the above expression can 
be integrated explicitly for arbitrary values of q. However, we find that the 
system admits oscillatory solutions only when q is a positive integer and for 
this choice we find the above expression reduces to either of the following 



forms depending on the value of q [15 



(i) Case 1 - q odd : We take q = 2m + 1, m — 0, 1, 2, Then equation 

flBT)]) reduces to the form 

x N{t) = r T : ; : 7 — (36a) 

[In- EL J cos ^ (sin 2 - f, + Efo C° k sin^"" 2 ) f ,)] 
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where £, = u t t + 5 h I = 1, 2, . . . , N and C° k = (2 k+1 m(m - 1) . . . (m - 
k))/{{2m - l)(2m - 3) . . . (2m - 2k - 1)). Using Q2ZD and (I55a|) we get the 
remaining N — 1 variables as 

*k(t) = , (36b) 

where fc = l,2,...,iV — 1. Note that for the solution (I36|) to be nonsingular 
periodic, we require the condition (71) / {wi)\Y^k=i ^1 + 1 < -^v- 

We note here that for the special case m = the solutions f )36ap and 
( 136bj) become the periodic/quasiperiodic solution 



x N (t) 



sin (t 



!n - Elli J cos (*,) 



4 sin (t 



**(*) = 7f~ V k = l,2,...,N-l. (37) 



which exactly matches with the solution given in for the so called coupled 
modified Emden equation. 

(ii) Case 2 - q even, Here we take q = 2m, m = 1,2,.... Then equation 
(I3"5|) reduces to the form 

x N (t) = — , (38a) 

[I N -D(t)+Bt]w 

where D(t) = Y^Li J cos^sin 2 ™" 1 U + Yl^i C e k sin^ 2 *-^ t,), 5 = (( 2 ™- 
l)!!)/(2 m - 1 (m - 1)!),' C L k = ((2m - l)(2m - 3) . . . {2m - 2k - l))/(2 fe (m - 
l)(m — 2) . . . (m — k)). Using (1271) and (I38ap we get the remaining N — 1 
variables as 

u\ 4 sin 4 
x k (t) = j— . (38b) 

[I N - D(t) + Bt]^+i 

From ( I38al) and (138bj) we find that equation (I30p admits only oscillatory 
dissipative type solution for the choice q = 2m, m = 1,2, . . . . 

Note that for q odd positive integer in ( 130]) . one can have either periodic or 
quasiperiodic solutions, depending on whether the uncoupled frequencies aVs 
are commensurate or not. In Fig. 1 we have presented quasiperiodic and peri- 
odic solutions for suitable choices of the uncoupled frequencies uj\, u 2 , ■ ■ ■ , ujn 
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with iV = 10 in the form of projected phase space plots in the X\ — £3 plane. 
We find that for the choice u\ = 2 and u 2 = CO3 = ■ ■ ■ = co w = 1, the system 
fl30|) exhibits 2 : 1 period oscillations. Similarly for the choice uj\ = \/2 and 
U2 = uj 3 = . . . = u w = 1, the system ( 130]) is found to exhibit quasiperiodic 
oscillations. 

One can note that equation (1501 can also be rewritten in the following 
first order form as was done by Iacono and Russo [lo| for the scalar case as 

Xi = UiUi - f(x)xi, (39a) 
iji = -u l x i - f{x)y u i = 1,2, . . . , N. (39b) 

Multiplying Eq. f)39ap by and Eq. f ]39bj) by x,- t and subtracting the 
resulting equations we get 

XiVi-XiVi = Vi(Xi +Vi). 
Dividing throughout by yf and rewriting we get 



1+ - , i = l,2,...,iV. (40) 



dt\UiyiJ \yi 

Upon introducing the angle variable 0, = tan _1 (xj/?/j), equation ()40l becomes 

ei = -u>i, t = l,2,...,N. (41) 

From (T4T|) it is obvious that the angle and hence the frequency (which is 
similar to that of the linear harmonic oscillator) are independent of the am- 
plitude of oscillation, irrespective of the form of f{x). However, one may note 
that this does not always imply isochronocity as the amplitude of oscillation 
may decay with time, as shown in equation (13"5|) . 

3. 2N integrals of motions 

In this section we show the existence of 2N independent integrals of mo- 
tion for equation (fl9"|) with / = /j, i = 1,2, . . . , N, being a homogeneous 
polynomial fT2"9l . In order to show that there exists iV time dependent inte- 
grals, let us consider the equivalent form of fT2"2"|) . that is, 

y A ^ ( Xi + f(x)xi) =cot ( Uit + Si ^ z = i ;2 ,...,iV. (42) 

Xi UJjXi 
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Figure 3: Projected phase space of system ([3"0]) in the x\ — X3 plane with N = 10, for two 
different values q — 3 (Figures (a,c)) and q = 5 (Figures (b,d)), respectively. Figures (a) 
and (b) describe the 2 : 1 period oscillations for the choice uj\ = 2 and 1^,^)3, ... , wio = 1. 
Figures (c) and (d) describe the quasiperiodic oscillations for the choice cji = v2 and 

W2,W 3 , ...,UW = 1. 



Rewriting the above equation (the first and the last expressions) we get N 
time dependent integrals as 



5; = cot 1 



Uit, i = 1,2,...,N, 



(43) 



where ^ = (fa + f(x)xi))/(ui). 

Now to find the remaining N integrals for q odd or even integer in (12TJ1) 
we proceed as follows. From (Ji2j) we get — = tan(6 l i ), 0{ = u{t + Si. Using 

this expression in the well known trignometric identity sin 9i = , tan 6i = and 

\J l+tan 2 9i 

cose,- cotc 



y/ 1+cot 2 6i 



we get 



sin(a;jt + Si) 
cos( + Si) 



X i 



X, 



Vi 



z = l,2, 



, N. 



x . 



x k sm(w N t+5 N ) 
xjv sin(u)it+ 

x\ sin 2 (wArt + S N ) 
x 2 N sin 2 (ojit + Si) 



Rewriting now (1271) in the form h = IkE^eJlI^I we obtain 



(44a) 
(44b) 

(45) 
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Now by using (I44ap we can write 



sin 2 ( c^t + Si 



2 i 2 ' 

Vi + xj 



sm 2 (u N t + 5 N ) 



UN + X iV 



(46) 



Substituting the above in (|4"5]) we obtain the first N — 1 time independent 
integrals as 



(4 + yl) o 



iV — 1. 



(47) 



For q odd positive integer, substituting ( 144")) and ( 147|) into ( I36aj) and 
rearranging we arrive at the following form for the integral 1^, 



T° 



1 



N 



{xn + Vn)*^ ' v/R+y?) V(x 2 + y 2 ) 

m— 1 



07 



yi 



X 



2m- 1 



2m-2fc-l 



fe=l 



(xf + yf) 



2m-2fc — 1 I J 



(48) 



which is the iVth time independent integral. For m = the integral (148 
becomes 



0*& + y&)* 



7"° 
'a 



(49) 



One can note that for m = the integrals (j43p . ( 14 7p and (149 p exactly match 
with those presented in For instance, equation (14"8"|) exactly reduces to 
the corresponding form given in 10] for the scalar case (N = 1). 

Similarly, for the case where q is even positive integer, we substitute ( |44"|) 
and (|47p into (I38aj) and rearranging we arrive at the following form for the 
integral Iff, 



1 



07 



X 



2m- 1 



m—l 



2m-2fc-l 



k=l 



A + yf, 



2m — 2fc — 1 



(50) 



We note here that the above first integral is a time dependent one. 
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4. Almost integrable systems 

In Sec. |3]we obtained the general solution of (|T9l) for the choice fi = f, 
i = 1, 2, . . . , N and / is a homogeneous polynomial as in (|29|) . However, we 
wish to point out that the system (1191) for any arbitrary choice of (1231) is 
almost integrable as there always exist N time dependent integrals (T4"3"j) and 
(N — 1) time independent integrals (j47|) . For complete integrability only the 
first order nonautonomous differential equation fT28|) needs to be integrated. 
For those forms of / for which this cannot be done explicitly, one can always 
carry out a numerical integration of fT28|) or apply a suitable approximation 
method to find Xjv(t). The remaining Xj(t)'s, i = 1,2, ...,(N — 1), can 
be obtained readily using f|2T|) and xjv(t). Equation fj4T]) ensures that if 
the solutions are periodic or quasiperiodic, the frequency is independent of 
amplitude. 




-0.04 -0.02 0.02 0.04 -0.2 -0.1 0.1 0.2 



X\ 

Figure 4: Projected phase space of the almost integrable system ([ST]) in the x\ — X2 plane 
for the choices (a) u>i = 1, u)2 = 2 exhibiting 1 : 2 period oscillation, (b) lo\ — \/2,u}2 = 1 
exhibiting quasiperiodic oscillation 

In order to demonstrate the above, let us consider the special case N = 2 
in ( 1301) and f(xi,X2) — x± + x\Xi + x\. This choice reduces equation ( 1301) to 
the following system of coupled second order nonlinear ODEs, 

i\ + [x\ + "ix\X 2 ^)X2 + (3xi + 3xiX 2 + 2x^X1 + x\x\(xi + 2X2 + 2^2) 

+2x\x 2 + xix Q 2 + x\ + c^xi = 0, (51a) 

X2 + (3xi + 3x1X2 + 2x 2 )xi + (x x + 3xix 2 )x2 + x x X2(xi + 2x2 + 2x 2 ) 

+xix 2 + 2x^x 2 + x^ + u\x 2 = 0. (51b) 
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Table 1: A comparison of the solution of Eq. (15^1) obtained using the variable step size 
Runge-Kutta fourth order (RK4) numerical procedure and the homotopy perturbation 
method (HPM). 



t 


x(t) using variable 


x(t) using HPM 




step size RK4 


up to third order 





0.2644 


0.2644 


1 


0.1556 


0.1884 


2 


-0.1416 


-0.1791 


3 


-0.2378 


-0.2647 


4 


0.1090 


0.1151 


5 


0.2381 


0.2846 


6 


-0.0128 


-0.0165 


7 


-0.2486 


-0.2932 


8 


-0.0790 


-0.0797 


9 


0.2663 


0.2824 



The solution of the above system of equations can be deduced by solving 
the following first order ODE obtained by the procedure discussed in the 
previous section, 

sm^cutf + di) sin [oJib + d{) 

However, we find that the explicit general solution of the above first order 
ODE is not known [14J. One can apply suitable numerical methods to solve 
this equation. FigureHJis plotted by solving Eq. ( 152|) using a variable step size 
Runge-Kutta fourth order method. We find the system (15T|) exhibits periodic 
and quasiperiodic oscillations which is shown in FiglH Here the projection of 
the phase space of (!5T|) in the X\ — x 2 plane for the set (u>i, u 2 ) = (1,2) and 
(tUijtc^) = (y/2,1) are shown for 1 : 2 periodic behaviour and quasiperiodic 
behaviour, respectively. 

One can also apply perturbation techniques such as the homotopy analysis 



16| to find approximate solution of high accuracy for Eq. f[52|) and compare 
the results with the numerical analysis. The homotopy perturbation method 
(l6i Il7l ] involves the introduction of an artificial parameter, say p, into the 
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original nonlinear equation A(u) = as 

(1 -p)(L(u) - L(u ))+pA(u) = 0, (53) 

where L is the linear operator corresponding to the linear part of the given 
equation and the A is the operator corresponding to the given nonlinear 
equation. Here uo is the lowest order approximate solution. Expressing the 
solution u as a power series in p, i.e. 

u = u + pui + p 2 u 2 + p 3 u 3 + . . . , (54) 

where «2,---, are the higher order approximations, one can substitute 
this series solution in (1521 to obtain a system of linear first order ordinary 
differential equations. Solving this system with suitable initial conditions 
one can obtain the approximate solution to the given equation in the limit 



p — > 1. For further details on this procedure one can refer to Refs. [16l . [17 . 

In Table [1] we compare the solution of Eq. (152]) obtained through this 
procedure with the numerical solution for the parametric choice 8i — 1, 82 — 
1.5, coi = y/2, uo 2 = 1. The values listed in Table Q] for the homotopy per- 
turbation method are calculated up to third order approximation. We note 
here that the accuracy of the perturbation solution will improve if further 
higher order approximations are taken into the calculation 17| , which we do 
not pursue here. 



5. Conclusion 

In this paper, we have shown that a system of N coupled nonlinear Lien- 
ard type oscillators admits periodic and quasiperiodic solutions or damped 
oscillatory periodic solutions with amplitude independent frequency of oscil- 
lations. We have derived explicit general solution and 2N integrals of motion 
for this system. Thus we prove this system to be completely integrable. We 
have also shown that another system of N coupled Lienard type oscillators 
is almost integrable in the sense that it admits (2N — 1) independent inte- 
grals and reduces to a single nonautonomous first order nonlinear differential 
equation. We have also shown that this almost integrable system also ex- 
hibits periodic and quasiperiodic oscillations for suitable parametric choices. 
For the general system of N coupled nonlinear oscillators f fl9|) with arbitrary 
form of nonlinearity, the nonlocal transformations reduce it to a system of 
iV coupled first order ODEs (1231) . It will be interesting to investigate further 
whether other forms of nonlinearity (different from fl29]) ) are also amenable 
to analysis. 
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